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Ivanov pointed out substantial analytical difficulties associated with self-gravitating, static, 
isotropic fluid spheres when pressure explicitly depends on matter density. Simplifications achieved 
with the introduction of electric charge were noticed as well. We deal with self-gravitating, charged, 
anisotropic fluids and get even more flexibility in solving the Einstein-Maxwell equations. In order 
to discuss analytical solutions we extend Krori and Barua's method to include pressure anisotropy 
and linear or non-linear equations of state. The field equations are reduced to a system of three 
algebraic equations for the anisotropic pressures as well as matter and electrostatic energy densities. 
Attention is paid to compact sources characterized by positive matter density and positive radial 
pressure. Arising solutions satisfy the energy conditions of general relativity. Spheres with vanishing 
net charge contain fluid elements with unbounded proper charge density located at the fluid-vacuum 
interface. Notably the electric force acting on these fluid elements is finite, although the acting elec- 
tric field is zero. Net charges can be huge (10^^ C) and maximum electric field intensities are very 
large (lO'^^ — IQ^"^ statvolt/cm) even in the case of zero net charge. Inward-directed fluid forces 
caused by pressure anisotropy may allow equilibrium configurations with larger net charges and 
electric field intensities than those found in studies of charged isotropic fluids. Links of these results 
with charged strange quark stars as well as models of dark matter including massive charged parti- 
cles are highlighted. The van der Waals equation of state leading to matter densities constrained by 
cubic polynomial equations is briefly considered. The fundamental question of stability is left open. 

PACS numbers: 04.40. Nr, 04.40. Dg, 04.20.Jb 



I. INTRODUCTION 

Self-gravitating fluid models are essential to many applications of general relativity, ranging from the historically 
important classical models of elementary particles to current computational and observational problems in astrophysics 
and cosmology. These systems are usually characterized by sets of physical variables which outnumber the independent 
field equations. Consistent problems can be posed if additional restrictions on the variables are specified, which may 
take the form of equations of state (EOS). 

Using an EOS to describe a self-gravitating fluid has important consequences when it comes to solving the field 
equations. For example, Ivanov [1\ has pointed out that finding analytical solutions in the static, spherically symmetric, 
uncharged case of a perfect fiuid with linear EOS is an extremely difficult problem. 

Interestingly, analytical difficulties may alleviate with increasing physical complexity. This situation has been 
illustrated by Sharma and Maharaj [2] in the case of a static, spherically symmetric, uncharged anisotropic fiuid. 
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These authors chose a particular mass function to reduce and easily solve the system of field equations combined with 
a linear EOS. 

The surprising simplification of the field equations for a charged perfect fluid satisfying a linear EOS was discussed 
by Ivanov [ij, who showed how to reduce the system involving the most general linear EOS to a linear differential 
equation for one metric component. However Ivanov also pointed out that the use of a polytropic EOS leads to 
non-integrable equations. 

Electrically charged fluids with anisotropic pressures constitute the next level of physical complexity. 

Charged, self-gravitating anisotropic fluid spheres have been investigated in general relativity since the pioneering 
work of Bonnor Q. Recently, this type of charged matter has been considered by Horvat, Ilijic and Marunovic in 
studies of gravastars [!]■ Models with prescribed EOS remain relatively unexplored. 

Motivated by MIT bag models of strange stars, Thirukkanesh and Maharaj [H derived solutions for charged 
anisotropic fluids with linear EOS from specific choices of one metric function and the electrostatic energy density. 
Their method - an extension of the procedure presented in for uncharged anisotropic fluids - provides motivation 
for the completely different approach developed in this paper. 

Current dark matter and dark energy models are associated with nonlinear EOS. Density perturbations may lead 
to fluid nucleation characterized by heterogeneous matter densities and pressures. Assuming specific EOS, bounded 
static fluid distributions and asymptotically flat spacetimes, we aim to study the arising self-gravitating objects. 
We insist on analytical solutions and note that prescribed EOS may cause serious trouble. However, we find that 
the convenient choice of a "frozen" internal metric reduces the Einstein-charged fluid equations to a system of linear 
algebraic equations for matter and electrostatic energy densities as well as anisotropic pressures. In a way, we are led to 
the simplest solution method available for this type of source with prescribed EOS and asymptotically flat spacetime. 
Models based on linear and nonlinear EOS are completely solved following essentially the same procedure. Remarkably, 
our approach uncovers the effects of different EOS on hydrostatic and electrical variables without interfering changes 
of internal metric (apart from adjustable numerical factors depending on junction conditions). This treatment offers 
a fresh view of the relationships among EOS, charge distributions and pressure anisotropy. 

In Section II we write the field equations and briefly review the work of Thirukkanesh and Maharaj. Also, the 
original Krori and Barua solution method is generalized to deal with charged anisotropic sources with regular interiors. 
Our approach to linear and nonlinear EOS leading to models with positive definite matter density is presented in 
Section HI. In Section IV we discuss junction conditions and reformulate our procedure in terms of dimensionless 
quantities. Section V includes detailed analysis of models with positive matter density and positive radial pressure, 
paying special attention to conditions for physical acceptability as well as equilibrium conditions and the values of 
key physical parameters in Gaussian-cgs units. In Section VI we speculate on the origin of charge in models with 
linear and nonlinear EOS, and check our solution method against a more complicated EOS describing quintessence 
stars. Finally we suggest avenues for further research involving stability analysis and gravitational collapse. 



II. CHARGED ANISOTROPIC MATTER 

The starting point is the static, spherically symmetric line element represented in curvature coordinates. It reads 

ds^ = e^df - e^dr^ - r^d0^ - sin^ 0d(l)^, (I) 

where v = ^{r) and A = A(r). For the static, charged source with density p — p{r), radial pressure Pr = Pr{r), 
tangential pressure pt — Ptir), proper charge density a — (j{r) and electric field E = E{r) the Einstein-Maxwell (EM) 
equations take the form 

87rp + i?2=e-^(^-l)+l, (2) 
8^^-i?2=e-^(^^ + l)-l, (3) 
8np, + E^^ — U + - + — (4) 
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where the primes denote differentiation with respect to r, and geometrized units (G = c = 1) are employed. 

Equations (H])-® are invariant under the transformation E — >■ —E, a — >■ —a. In this work we exclusively deal 
with the positive square root of E^ . 

With our choice of radial coordinate the metric function and electrostatic energy density E"^ assumed in [H take 
the forms 

1 + (a - 6) r2 ' 



E = -7^. (7) 



(1 + ar2) 



where a, b, k are arbitrary constants. 

The combination of ^ and ([7]) with ^ yields 



87r(l + ar2)' 

Hence p and E"^ are proportional when k ^ and k ^ b. This expression for p(r) is joined to a linear EOS to 
provide Pr{r). The explicit form of together with ([6]), (O and ([3]) imply a linear differential equation for iy{r) 

which is analytically solved. Finally, when ^ and the explicit form of ^{r) are substituted into ([4]) we get the 
corresponding expression for pt(r). 

The above procedure leads to analytical solutions which depend on a number of free parameters. The analysis 
presented in Q considers sensible choices for these parameters in the context of charged stars. However we find a 
singularity in the charge distribution at r = when Q and ^ are combined with ([5]). 

From ([6]) we see that A(r) satisfies A(0) = 0. Assuming k > 0, ^ implies 



— (9) 

for small r, which diverges at r = 0. The electric field associated with ([7]) does not vanish at r = 0. This choice for 
E{r) prevents the regularity of the charge distribution at the centre of the sphere. 

In our view, the vanishing of the electric field at the center of a spherically symmetric charge distribution should 
be a condition for physical relevance of the solution. Furthermore the significance of charge density singularities is 
unclear and we aim to investigate their possible occurrence in models satisfying E(0) = 0. To this end, we extend 
Krori and Barua's approach to charged isotropic fluid sources @ and deal with charged anisotropic sources with 
prescribed EOS. 

Krori and Barua (KB) constructed singularity-free models of static, charged perfect fluids with metric ([T]) given by 

A = Ar2, (10) 

p = Br^ + C, (11) 

where A, B and C are constants. This internal metric satisfies the conditions for regularity at r = discussed by 
Lake and Musgrave 0. As a consequence of this choice, equations ©-Q with Pr = Pt = P were reduced to a system 
of three linear algebraic equations for p, p and E'^. Furthermore a was obtained combining ^ with the chosen square 
root of E^ and the assumed form of A. 
We re-write the field equations as 

(12) 
(13) 
(14) 

where f{r),h{r),j{r) are determined by the right sides of (H))-©, (flO)) and dTTI) namely 

f{r)^e-^^' (2A-\)+\, (15) 



8np^ 




= fir), 


Snpr - 


-E' 


= Hr), 


Snpt - 


hE^ 
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hir) = e— - ( 2S + 1^ ^ 



j(r-) - e-^''' [S + B^r^ + {B - A) - ABr^] . 
We impose center and boundary conditions on E{r) and Prir) respectively: 

E{0) = 0, 

Pr{a) = 0, 

where a is a positive constant and r ~ a defines the charged fluid-vacuum interface. 
From (inil-dlll), (HSll-dlTl), ^ and ^ we obtain 

87r/9(0) = 3A, 



47rp(a) = {A + B)e 



87rpr(0) = 2B- A, 



(16) 
(17) 

(18) 
(19) 

(20) 

(21) 
(22) 
(23) 



87rpt(a) = e" 



4B - A + aMS^ 



AS +- 



(24) 



2B 



(25) 



From (22) and ([25]) we conclude that only one pressure value is associated with r — 0. Hence Pr{Q) and pt(0) denote 
the same quantity i.e. central pressure. 

Equations p ^ -p ^ yield general expressions for pt and E'^ namely 



j{r) - /(r) 
8. 



E^ = fir) - Snp, 



where p is still undetermined. 



(26) 



(27) 



III. LINEAR OR NON-LINEAR EQUATION OF STATE 

At this point we select an EOS with the general form 

Pr = Pr(p,ai,a2), 

where ai and a2 are constant parameters. 

These parameters are constrained by the system of equations 



(28) 



Pr{0) =Pr [p{0),ai,a2] , 



(29) 



= [p{a),ai,a2] 



(30) 
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Combining and ([T3)) we get 



fir) + /^(r) 

P + = g:;;^ ' (31) 



which may be solved with the assumed EOS to generate specific forms of p and pr. The arising p is put into (j26p and 
([?f)) to yield general expressions for pt and i?^. 
Firstly we consider the linear EOS 

Following the outlined procedure we obtain 



Pr ^ ai+ a2P- (32) 
[/(r) + hir)] ~ ai 



1 + a2 



(33) 



1 + Q!2 

j(r) + h{r) + aa [j(r) - /(r)] - 87rai 

ft = TTTT — A ' 35) 

87r(l + 02) 



1 + ^2 ■ ^ ^ 



Taking into account ([5^ we solve the system ([^ - (15(1)) and get 



- - p{0) - p{ay "^-p(0)-p(a)- 



Secondly we deal with the nonlinear EOS 

ft = /3i + ^, (38) 

where n 7^ — 1. Combining psp with (1311) we obtain a polynomial equation for p which is quadratic only for n — —2 
or n — 1. The second choice for n yields 

fc(r) - /3i ± 7[A;(r)-/?i]'-4/32 
P= ^^^^ , (39) 



where 



M.) ^ ^fc^. (40) 



Solving dMl) and ^ with ([381) and n = 1 we get 



p(0)p.(0) p(0)p(a)p.(0) 

p(0)-p(a)' p(0)-p(a) ■ ^ > 

li P2 < then each root in p9p has definite sign. Particularly, the root with the positive radical determines a positive 
definite matter density. 

Equation (pS)) is a modification of the Chaplygin gas EOS used by Bertolami and Paramos to describe neutral 
dark stars [8]. Actually, the Chaplygin gas EOS is polytropic with negative polytropic index. The additional term /?i 
allows for a charged fluid- vacuum interface where Pr = 0- 

Thirdly we consider the modified Chaplygin EOS 

Pr = 7iP + — (42) 
P 
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used in the study of static, neutral, phantom-like sources presented by Jamil, Farooq and Rashid @. 
Equations (HH) and combined with (gl]) yield 



7i 



p(0)p,.(0) 



p(0)2-p(a)2' p(0)2-p(a)2- 
Putting (|42]) into ((3T|) we get a quadratic equation for p which admits the solutions 



p{0)pr{0)p{ay 



fc(r)±Vfc(r)2^ 4 (1 + 71)72 
2(l+7i) 



(43) 



(44) 



If 1 + 7i > and 72 < then the sign of the radical in fj44l) uniquely determines the sign of p. As in the previous 
case, a positive radical leads to positive definite matter density. 

Finally, the expressions for p selected from (15^ and (H^ are combined with the corresponding EOS as well as ([25)1 
and (|27)) to yield solutions for pr, pt and E"^. 

It is clear that the constant parameters appearing in psp play a central role in our extension of the KB method. 
Their values determine the existence of positive and negative definite matter densities p{r). Furthermore they lead to 
numerical values for quantities like , which is associated with sound propagation in anisotropic charged fluids [lo| . 



IV. JUNCTION CONDITIONS AND ADIMENSIONAL FORMULATION 



Equations (fBT)) . (|4T|) . (j43| show that the constant parameters are determined by p{0), p{a) and Pr(0), which, 
according to (PH)) and (P^ . are functions of the KB constants A and B appearing in pUj) and (ITT|) . Notably, these 
two constants as well as C are fixed by suitable junction conditions imposed on the internal and external metrics at 
r = a. 

The external Reissner-Nordstrdm (RN) metric is given by ([T]) with 



= 1 - 



2m 



,A(r) 



2m q 



2 1 



(45) 



Junevicus [T2[ derived expressions for A, B, C from the continuity of the first and second fundamental forms across 

the surface of the charged fluid sphere. In terms of the dimensionless parameters = — and x = ~ tiis results take 
the form 



A = 



In (1 - 2/i + 



(46) 



B = — 



M - X 



1 - 2m + 



(47) 



C 



In (1 - 2/i + X 



(48) 



1 - 2/i + 

Combining ((25|) with (|46p and (|47)) we obtain the electrostatic energy density at the surface of the sphere namely 

„2 



E{aY 



(49) 



On the other hand, using pS)) in the electrovac {p = 0) case of ^ we see that the same quantity evaluated in the 
external (r > 0) spacetime region is 



(50) 



These results are compatible with the continuity of the electric field at r = a. 

For simplicity in numerical calculations, we reformulate our results in terms of dimensionless quantities. From (jlOp 
and (fTTI) we see that A^ B have dimension of length~^ whereas C is dimensionless. Clearly /(r), h{r), j{r) as well 
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as p, pr, pt, E"^ , (J also have dimension of length~'^. We get the dimcnsionlcss version of any variable or constant 
parameter multiplying by the appropriate power of a. Here adimensionality is denoted by tildes, though quantities 
which arc originally dimensionless are denoted by the original symbol. The dimensionless radial marker a; = ^ is used 
so the interior of the fluid sphere is described with x & [0, 1). 
For the first model we get 



Pr = 011+ 0L2P, 



(51) 



1 



/>) + K^) 



1 + a2 



(52) 



Pr 



1 + Q!2 



(53) 



where 



Pt 



j(.t) + //(.r) +02 %r)-f(.r) - Sto 



87r(l + 02) 



-2 _ STTcii + a2.f{x) - h{x) 



1 + a2 



p(l)p.(0) 
p(0) - p(l) 



Pr(0) 



p(0)-p(l)' 



(54) 

(55) 

(56) 

(57) 



ft(a;) = e-^=^ \2B 



1 \ 1 



(58) 



2B + 2B'^x^ + 2 



(B-i) 



2ABx' 



m = 



3A 
8^' 



(59) 
(60) 



4-77 



(61) 



The second model is described with 



Pr(0) = 



2B-A 



A = a^yl, B = a^B. 



Pr= Pl + —, 
P 



(62) 
(63) 

(64) 
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kix)-(3i±d -4/3: 



~ mPAO) r. p(0)p(l)p.(0) 
"l = "^77T^ ^TTT' P2 - 



For the third model we obtain 



Pr = 7lP+ — , 
P 



(65) 
(66) 

(67) 



fc(x) ± ^fc(x)2- 4 (1 + 71)72 



2(1 +71) 



(68) 



71 



P(0)Pr(0) 
p(0)2-p(l)2^ 



72 = 



p(0)p.(0)p(l)2 

p(0)2-p(l)2 



(69) 



For models with nonlinear EOS, pt and i?^ are calculated with the dimensionless versions of ((26)) and (|27)) respec- 
tively. 

From (j49|) we get 



(70) 



which applies to the three models. 

Finally the dimensionless proper charge density is given by 

6 2 0! 

Attx-^ dx 



(x^e) 



(71) 



V. ANALYSIS OF THREE MODELS 



We consider only external solutions (|45|) which exclude horizons. The standard analysis of the roots of 500 — 
implies that the values of x ^re restricted by the values of /x. Three cases arise: 

1. A*e(0,i), 2^-i<x'<m'; 

2. ^e(0,i), x^p] 

3. > 0, X > P- 

From equations and we see that p{0) is positive if and only ii 2^ — 1 < x^ < "^p, which we admit as an 
additional restriction on x- 

The geometric mass in (H5t is given by m = where G is the gravitational constant, M is the mass of the 

charged source, and c is the speed of light - all in conventional units. If the mass M of the charged star equals one 
solar mass and a = 10 Km, then /i = w 0.147. Hence 2// — 1 is negative and y/2fl « 0.543. 

Using (l60l) and (|62|) the central density p(G) and central pressure Pr{0) can be evaluated for /i — 0.147 and 
X G [0, 0.543). The monotonic decreases of these two parameters with increasing x are shown in Figures 1 and 2. We 
observe that Pr{0) is positive only when x takes values in the range < x < xoj where x = Xo is the only zero of 
Pr{0) in the interval x € [0, 0.543). It is approximately given by xo ~ 0.191. 

For fixed we see that p(0), p(l) and Pr{0) depend only on x- Hence equations ([56| . (|66t and (|69)) provide 
expressions for the EOS parameters as functions of x- As shown in Figures 3-5, $2{x) < 0, 71 (x) > and 72 (x) < 
for /i — 0.147 and x G [0,Xo)- These results allow us to select the positive definite roots from (|39|) and ([44]) . and 
discuss models with p(a;) > 0. 

For concreteness we carry out the numerical and graphical analysis of models with /i = 0.147 and x € [OjXo); 
characterized by positive matter density and positive pressure at a; = 0. Clearly x = describes sources with zero net 
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charge. This condition is compatible with non-zero proper charge densities a{x). Higher values of x ^^re associated 
with increasingly repulsive electrostatic forces that affect pressure and matter density profiles. We aim to find out 
how the different forces which allow equilibrium configurations accommodate to varying net charge. 

To begin with we examine sources with linear EOS and a selection of x values. The corresponding profiles for 
matter density, radial and tangential pressures, and electric field are displayed in Figures 6-9. We observe decreasing 
values of matter density and pressure associated with increasing values of x- Figure 9 shows electric field profiles 
satisfying £'(0) = and E{1) = X- Very large negative gradients of E develop near the surface of the sphere for 
X = 0, and E increases at every x € (0, 1] with increasing x- Notably Figure 7 indicates that increasing net charges 
and electric fields are associated with decreasing radial pressure gradients at each point. 

Figures 10 and 11 show profiles for the squares of radial and tangential sound speeds, defined by v"^^ = ^ and 

^s* = ^ respectively. We observe that v^^. is independent of x and decreases with increasing x- For fixed Xj ^st 
monotonically decreases with increasing x; and increases with increasing x ioi fixed x. These parameters satisfy the 
inequalities < < 1 and < < 1 everywhere within the charged fluid for the five values of x considered. 

Based on the standard analysis of energy conditions for charged anisotropic fluids (see, for example, ^Ifl]), we find 
that the weak energy condition (WEC), the strong energy condition (SEC) and the dominant energy condition (DEC) 
are simultaneously satisfied if and only if the following six inequalities hold at every point within the source: 

P+Pr>0, (72) 

P+^>0, (73) 

p-i-pt-f— >0, (74) 
47r 

E? 

p-|-p^-|-2pt + — >0, (75) 
47r 

EP' & 

/^+^-|p. + ^l>0. (77) 

Inequalities (|72|) and ((73|) hold automatically for the sources considered here. Straightforward plotting of the left 
sides of ((74l) - ((77)) shows that these inequalities are satisfled as well at every x e [0, 1]. 

Further analysis shows that models with nonlinear EOS (|64|) and (|67| are very similar to those satisfying ()51|) . 
Particularly, matter densities and radial pressures are everywhere positive, matter densities as well as radial and 
tangential pressures monotonically decrease with increasing x, radial and tangential sound speeds satisfy < < 1 
and < < 1, VsT decreases with increasing x for fixed x, and WEC, SEC and DEC are satisfled for /i = 0.147 and 
X = 0,0.05,0.1,0.15,0.19. However models with nonlinear EOS get Vsr proflles which increase monotonically with 
increasing x and fixed x- 

Another difference among the three models concerns the anisotropy parameter A = — p^.. As shown in Figures 
12 - 14, each EOS affect the dependence of A on x in a different way. Particularly, the sign of A is the same at every 
X e (0, 1] for fixed x only in models with linear EOS. The effect of on pressure anisotropy is notorious as three 
A(a;) profiles develop sign changes in that case. However the predominance of radial pressure over tangential pressure 
for the highest value of x is a common feature of the three models. 

Are these three types of models physically meaningful? Delgaty and Lake found that only a relatively small number 
of proposed perfect fluid sources for the Schwarzschild metric satisfy a set of well established conditions for physical 
acceptability [13]. These conditions include regularity of the origin, positive matter density and pressure, decreasing 
matter density and pressure with increasing r, causal sound propagation and smooth matching of the internal and 
external metrics at the fluid- vacuum interface (r = a). Some authors add the condition of monotonically decreasing 
sound speed with increasing r (see, for example, [l^). Regarding anisotropic fluid sources, the causality condition 
< < 1 has been imposed on radial and tangential sound speeds (see, for example, [lH| and [3). Disagreement 
arises in connection with the sign of tangential pressure which is unrestricted for some authors (see, for example. 



10 



[isj ) and strictly positive for others (see, for example, [iBl)- The above discussion indicates that our three charged 
anisotropic sources satisfy most of the usual acceptability criteria. Conflict may arise only in connection with negative 
tangential pressures occurring for the highest x values, constant sound speeds Vsr associated with the linear EOS, 
and the increase of v^r with increasing x appearing in the models with nonlinear EOS. 

Electric interactions and charge distributions in our models deserve further analysis. We have found that E profiles 
are affected by the choice of EOS. Particularly, in models with x = the maximum E values are approximately 0.03, 
0.07 and 0.05, corresponding to (|5ip . ([M)) and (|57)) . respectively. In the case of EOS ([5T|) . (t(0) is finite and increases 
with increasing x; and (t(1) is finite for x — 0.05, 0.l, 0.15, 0.19 but unbounded for % = 0. Moreover, only sources with 
X = contain electric charge of both signs. We could have anticipated this behaviour of a from Maxwell equation ((7T|) 
and the slope changes displayed in Figure 9. Figure 15 shows the arising a profiles in the restricted interval [0, 0.999]. 
The most important difference among a profiles associated with the three EOS is that models with nonlinear EOS 
and X — 0-05 also involve positive and negative electric charges. 

We have noticed that the charge distributions considered by Thirukkanesh and Maharaj (H are singular at the 
origin, where the electric field does not vanish. Alternatively, all our sources have vanishing electric field and finite 
proper charge density at r = 0. On the other hand, our models with x = involve charge distributions which 
are singular at the fluid- vacuum interface. We remark that charged thin shells are not considered in our study and 
electric fields are continuous at r = a. Also, charged sources with x = are characterized by E{a) = and unbounded 
limr_^£t- E'{r). Arbitrarily large electric field gradients near the fiuid- vacuum arc puzzling and we proceed with a 
preliminary discussion of their significance. 

The proper charge density a appears explicitly in the inhomogeneous Maxwell equation ([5]). It determines the net 
charge inside a sphere of radius r through the formula 

f A(t-) „ 

q{r) = An a{r)e = r dr. (78) 

^0 

We remark that q{a) is the total charge of the source, denoted by q in (|45|) . It is q{r) the quantity that actually 
determines the electric field 

E{r) = ^. (79) 

Proper charge density (j{r) is not an observable quantity in standard EM theory. We have seen that E{r) takes 
part in the formulation of the energy conditions of general relativity through the contribution of the Maxwell field 
to the energy-momentum tensor. In contrast, the authors are not aware of any condition for physical acceptability 
imposed directly on a. 

Both spacetime curvature and electric field gradients deviate the worldlines of charged test particles. Balkin, van 
Holten and Kerner [l3| and van Holten [l^ have derived covariant formulae for the relative acceleration of particles 
with the same charge/mass ratio. Their equations include terms containing derivatives of the Maxwell tensor. The 
question arises as to whether these terms could lead to infinite relative accelerations for pairs of charged test particles 
passing through the fluid- vacuum interface of models with x = 0. Hence very large relative accelerations of test 
particles could allow indirect observation of inflnite charge density located at the vacuum-fluid interface of our sources. 

The fundamental invariant of the electromagnetic field Ii = F^F^^'^ has been used in the analysis of genuine 
singularities of static solutions for the Einstein-Maxwell equations [191 . This invariant is bounded at the fiuid- vacuum 
interfaces of our three types of models. The inhomogeneous Maxwell equation implies that I2 — F^'^ ■i.F'^ ^^.^ is 
proportional to o'{r)'^- The fact that <j{r) is unbounded at r = a in models with x = determines the singular 
behaviour of a differential invariant of the theory. Questions about physical acceptability of solutions with regular Ji 
and singular I2 are analogous to those regar ding spacetimes with regular curvature invariants and singular differential 
invariants discussed by Musgrave and Lake |2q |. 

The discussion of physical acceptability for sources with zero net charge points at the equilibrium condition for 
sections with infinite charge density. Fluid elements with unbounded cr are located at r = a, where the electric field 
vanishes due to the choice x = 0. Nevertheless, the conclusion that no electric force acts on elements of charged 
fluid with infinite a is not straightforward and deserves closer examination. Other features of our models motivate 
further analysis of equilibrium conditions. For example, as we keep /i constant and increase Xj repulsive electric 
forces increase with decreasing matter density, decreasing pressure and decreasing radial pressure gradients. The 
question is, how gravitational and other fluid forces counteract increasing electrostatic repulsion when the charged 
fluid becomes more diluted and pressure gradients tend to vanish? The situation is clarified using the generalized 
Tolman-Oppenheimer-Volkov (TOV) equation as presented by Ponce de Leon [l^]. (See also [HI.) It reads 



Mn (p + Pr) h:z±. dpr OA 2 

dr r 



- ^ + ^^e* + r - Pr) = 0, (80) 
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TABLE I: Dimensionless values of physical variables for the model with linear EOS. 



X 


Po 


Pa 


Pta 


E-i 







0.042 


0.003 


0.001 





0.0009 


0.19 


0.036 


0.00005 


-0.002 


0.036 





where Mq = M^^r) is the effective gravitational mass inside a sphere of radius r and q = q{r) is given by ((75)) . The 
effective gravitational mass is given by the expression 

Mair) = Ir^e^iy', (81) 

derived from the Tolman-Whittakcr formula and the Einstein-Maxwell equations [l^l- 

Equation (|80|) expresses the equilibrium condition for charged fluid elements subject to gravitational, hydrostatic 
and electric forces, plus another force due to pressure anisotropy. Combined with ([TOl) . (jlip . (|63|) and ([79]) . it takes 
the adimensional form 

Fg + F,, + Fe + Fa = 0, (82) 

where 

Fg = -Bx {p + pr) , (83) 

Fe^aEe^, (85) 
Fa^- {pt - Pr) ■ (86) 

X 

The profiles of Fg, Fh, F^ and Fa for sources with linear EOS, /i — 0.147 and x = are shown in Figure 16. Notably, 
the electric force acting on fluid elements with unbounded (t located at a; = 1 is finite although E{1) — in this case. 
Equivalently limj^^i aE is flnite. Furthermore F^ is the weakest force and changes sign at a; « 0.85. Both Fh and Fa 
point outwards at every x € (0, 1]. Electric forces increase and radial pressure gradients decrease with increasing x- 
When X = 0.19 the hydrostatic force F/i is relatively insignificant and Fa points inwards, acting in conjunction with 
gravitational attraction to compensate the electrostatic repulsion. This situation is described in Figure 17. Clearly, 
the sign of Fa changes due to the predominance of Pr over pt for the largest x values. This sign inversion is essential 
for the configuration of our static, charged anisotropic sources with linear EOS. 

We have extended the analysis of these four forces to models with nonlinear EOS and found essentially the same 
equilibrium configuration discussed above for linear EOS and x — 0.19. Particularly, Fh is comparatively small, Fg 
and Fa have nearly the same profile and jointly counteract the electrostatic repulsion. Hence the choice of EOS has a 
negligible effect on the compensation of relatively large electrostatic repulsion by gravitational attraction and pressure 
anisotropy. 

Equations ©-([S]) indicate that matter density, radial pressure, tangential pressure and electric field strength affect 
the spacetime metric in our relativistic fluid models. Our choice /x = 0.147 corresponds to a compact stellar object 
(M = 2 X 10'^^ g, a = 10^ cm), so pressure is expected to play a substantial role here. The question arises as to 
whether the electrostatic energy density signiflcantly contributes to the source of gravity. We compute dimensionless 
values of the physical variables and get maxima and minima for po (central matter density), pQ (central pressure), pta 
(tangential pressure at r = a) and F^ (electrostatic energy density at r = a) as well as F^ (maximum electrostatic 
energy density for x — 0)- Approximate numerical results for the model with linear EOS are displayed in Table I. 

We observe that the maximum values of po and pta are just one order of magnitude smaller than the maximum of pq. 
These maxima correspond to the source with x = i.e. zero net charge and weakest electrostatic repulsion. Notably, 
these density and pressure values are very similar to the dimensionless density ps and dimensionless central pressure 
Pos of the (uncharged) Schwarzschild internal solution (SIS) with the same /i value. Using well-known formulae for 
the SIS in geometrized units (2]| we derive the dimensionless expressions 
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TABLE II: Density and pressures for the model with linear EOS expressed in cgs units. 



X 


po {gr/cm^) 


po {dyn/cm'') 


Pta {dyn/cm'') 





5.6 X 10"'' 


3.6 X 10^" 


1.2 X 10^" 


0.19 


4.8 X 10'* 


6 X 10^^ 


-2.4 X 10^" 



TABLE III: Electric field strengths and net charge for the model with linear EOS expressed in Gaussian-cgs units 



X 


Ea (statvolt / cm) 


Em {statvolt / cm) 


q (statcoul) 








3.7 X 10^^ 





0.19 


2.3 X lO""* 




1.9 X 10^"* 



Pos^^, (88) 

leading to the approximate results ps = 0.036 and pos — 0.003. 

The value of E^^ shown in Table I indicates a substantial contribution of the Maxwell field to the source of gravity 
in the case of zero net charge, although matter density and pressure predominate. 

In sources with maximum net charge {x — 0.19) po decreases in two orders of magnitude, pta doubles its absolute 
value and gets the opposite sign, and po and have the same approximate value. From figures 13 and 14 we see that 
the reduction of radial pressure with increasing net charge makes gravitational attraction weaker in spite of the larger 
values of electrostatic energy density. We have discussed above that the sign inversion of tangential pressure leads to 
an extra force which complements gravitational attraction, so that the stronger electric repulsion gets compensated. 
Certainly static equilibrium is attainable in these charged sources thanks to pressure anisotropy. 

The discussion of the above quantities in conventional units is interesting as well. Adimensional density and pressure 
values are expressed in geometrized units and then converted to cgs units using well-known conversion factors [2^. 
Similarly, adimensional values of electrical variables are expressed in geometrized units, then converted to Heaviside- 
Lorentz units and finally to Gaussian-cgs units. The corresponding results are displayed in tables II and III. 

These values of central matter density and pressure as well as electric field strength and net charge in conventional 
units are similar to those discussed by previous authors in the context of charged compact objects with isotropic 
pressure [1^. It is understood that the huge gravitational attraction determined by these values of central matter 
density and pressure compensates repulsive electrostatic forces associated with field strengths and net charges with 
orders of magnitude 10^* statvolt / cm {10^^ V/ cm) and 10^^ statcoul (10^® C) respectively. However our study of 
charged anisotropic sources points at the crucial role that forces caused by pressure anisotropy can play in the 
construction of equilibrium states. This anisotropic effect is equally important for the three models considered here, 
at least for the assumed value of fi = 0.147. Provided that pressure anisotropy supplies an inward-directed force that 
compensates electrostatic repulsion, we expect our sources can achieve higher charges and electric field strengths than 
sources with isotropic pressures and the same p value. Also, our results for sources with zero net charge suggest a 
possible role for pressure anisotropy and electric charge in the structure of static sources for the external Schwarzschild 
metric. 



VI. CONCLUDING REMARKS 



The linear equation (|32|) links our analytical approach with the numerical treatment of electrically charged strange 
quark stars by Negreiros et al. [l^l based on the MIT Bag Model. They also assume (|18p and impose vanishing 
isotropic pressure at the fluid- vacuum interface. 

We have borrowed nonlinear EOS ([55]) (with n — 1) and (UU from current models of dark matter and dark 
energy. Most applications of Chaplygin and modified Chaplygin gases are cosmological and describe non-static, 
neutral gravitational fluids with isotropic pressures. Instead our bounded sources involve static, asymptotically flat 
spacetimes and charged anisotropic fluids. Actually we have targeted the effects of these EOS on the interior regions 
and fluid- vacuum interfaces of charged stars appart from any cosmological framework. The question arises as to how 
dark matter-dominated star interiors could get electric charge. Speciflc charge transfer mechanisms hypothesized in 
studies of charged neutron stars [25^ could be considered. Alternatively, a fraction of dark matter could be made 
up by massive particles with electric charge (CHAMPs) [26.], so dark stars could be charged. We remark that our 
procedures and results are totally independent of any charge generation mechanism. 

Lobo initiated the study of van der Waals (VDW) quintessence stars [131 . Apart from variations of p and Pr in 
the interior of the VDW fluid, which may be caused by gravitational instabilities, this author imposes a cut-off of the 
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energy-momentum tensor at r = a, where the internal metric matches the external Schwarzschild solution. Again, 
the motivation for introducing this new type of bounded source is cosmological. Provided that the quintessence EOS 
leads to interesting descriptions of the late universe, stellar objects arising from fluid nucleation through density 
perturbations are explored. 

The extended KB approach can be applied to a sphere of charged quintessence fluid. Apart from the hypothetical 
link of dark matter with CHAMPs, we cannot justify the introduction of charge in this model. However we aim to 
test the applicability of the present approach against the more complicated VDW EOS 

Pr^j^-5,p^ (89) 
1 - dzp 

which describes dark matter and dark energy as a single fluid [27l| . 

Assuming that the interior metric is joined to (j45l) . and that the electric field and radial pressure vanish at r = and 
r — a respectively, we evaluate (|89p at the center and boundary of the charged sphere and solve two linear equations 
for 5i and 82 as functions of p(0), p(a), Pr(0) and 5^. Then ([5^ is combined with ([5T|) to obtain 

|52(53P' - (^2 + &:^)p' + + 63 .g(r)] p - g{r) = 0, (90) 

where g(r) = ^^^±tMr). 

We do not go further into the analysis of the real roots of (|90l) and the corresponding sources. However, we anticipate 
difficulties in the selection of /z, x ^-iid '^3 values leading to positive definite matter densities. 

The analysis of stability is beyond the scope of this paper. Results presented by Andreasson regarding stability 
properties of charged anisotropic spheres could shed light on this fundamental issue. We highlight the potentially 
interesting study of gravitational collapse of sources with zero net charge, and hope our results will motivate further 
research on these topics. 

Finally we point out that different EOS could be matched at specific radii to provide composite models. For 
example, a linear interior (core) with a nonlinear exterior (layer) matched at a location inside the charged fiuid sphere 
might provide a better model than a simple EOS throughout. Using the extended KB approach developed here, 
both the core and the layer would be described with the same general metric ([T|) and the ansdtze (fTO|l and (jlip with 
presumably different values of A, B and C in each region. Less restricted choices of these constants could lead to 
charged thin shells emerging at the interface between the fluids. 
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FIG. 1: Central density p(0) as a function of x for M — 0.147 (y = p(0)). 




FIG. 2: Central pressure PriQ) as a function of x for M = 0.147 (y = PriQ))- 
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FIG. 5: Coefficient 72 as a function of x for M = 0.147 {y = 72). 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 6: Matter density p as a function of x for five different values of x- In this figure the curves 1 — 5 from the top corresponds 
to X = 0, 0.05, 0.10, 0.15 and 0.19 respectively {y = p). 
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FIG. 8: Tangential pressure pt as a function of x for five different values of x- In this figure the curves f — 5 from the top 
correspond to x = 0, 0.05, 0.10, 0.15 and 0.19 respectively {y = pt). 
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FIG. 10: Square of radial sound velocity v^r as a function of x for five different values of x- In this figure the curves 1 — 5 from 
the top correspond to x = 0, 0.05, 0.10, 0.15 and 0.19 respectively (y = Vg^). 
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FIG. 12: Anisotropic parameter A for the first model as a function of x for five different values of x- lu this figure the curves 
1 — 5 from the top correspond to x = 0, 0.05, 0.10, 0.15 and 0.19 respectively (y = A). 
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FIG. 14: Anisotropic parameter A for the third model as a function of x for five different values of x- In this figure the curves 
1 — 5 from the top correspond to x = 0, 0.05, 0.10, 0.15 and 0.19 respectively (j/ = A). 
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FIG. 16: Four different forces acting on fluid elements in static equilibrium as functions of x for /i = 0.147 and x — 0. In this 
figure the curves 1-4 from the bottom (blue, green, tan, cyan curves) correspond to Fg, F^, Fa, F^, respectively. The dummy 
variable y on the vertical axis represents any of these forces. 
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FIG. 17: Four different forces acting on fluid elements in static equilibrium as functions of x for = 0.147 and x = 0.19. In 
this figure the first (green) and second (cyan) curves from the top correspond to Fe and Fh, respectively. The third (tan) and 
fourth (blue) curves intersect at a; ~ 0.75 and represent Fa and Fg, respectively. The dummy variable y on the vertical axis 
represents any of these forces. 



